Magnetic correlations in YBaCo407 probed by single-crystal neutron scattering 
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We have studied the frustrated system YBaCo407.o generally described as an alternating stacking 
of Kagome and triangular layers of magnetic ions on a trigonal lattice, by single crystal neutron 
diffraction experiments above the Neel ordering transition. Experimental data reveals pronounced 
magnetic diffuse scattering, which is successfully modeled by direct Monte-Carlo simulations. Long- 
range magnetic correlations are found along the c-axis, due to the presence of corner-sharing bipyra- 
mids, creating quasi one-dimensional order at finite temperature. In contrast, in the Kagome layers 
(afe-plane), the spin-spin correlation function -displaying a short-range 120° configuration- decays 
rapidly as typically found in spin-liquids. YBaCo407 experimentally realizes a new class of two- 
dimensional frustrated systems where the strong out-of-plane coupling does not lift the in-plane 
degeneracy, but instead act as an external "field". 

PACS numbers: 75.30.Kz, 72.80.Ga, 25.40.Dn 
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Geometrically frustrated (GF) magnets continue to at- 
tract both experimental and theoretical attention due 
to the fact that exotic magnetic states such as spin-ices 
and spin-liquids emerge from frustrated magnetic inter- 
actions. In such systems, the impossibility of uniquely 
minimizing all energy terms in the magnetic Hamilto- 
nian, leads to a macroscopically degenerate manifold of 
ground states, which in the absence of perturbations, pre- 
vents the formation of long-range magnetic order. A clas- 
sic example of GF systems extensively studied theoreti- 
cally PI include vertex-sharing triangular or tetrahedral 
networks with antiferromagnetic (AFM) interactions, re- 
alized experimentally in, for instance, SrCr8Ga40ig 

The recently discovered class of compounds RBaCo407 
R=rare-earth, Y) has attracted considerable attention 
0] since it is among the few crystalline solids whose 
structural motif embraces magnetic Kagome layers, a 
classical GF lattice. The transition metal magnetic sub- 
lattice of RBaCo407 is made of triangular and Kagome 
layers of Co ions (Co(l) and Co(2), respectively) that 
are alternatively stacked along the c-axis of the trigonal 
crystal structure, as shown in Fig[TJ Strong AFM in- 
teractions are mediated through Co-O-Co superexchange 
pathways between first neighbor Co ions, tetrahedrally 
coordinated by oxygens, creating an exchange topol- 
ogy based on triangular motifs. Magnetic susceptibility 
measurements on all analogs of the RBaCo407 family 
are characteristic of strongly frustrated systems with a 
building up of strong correlations much above the or- 
dering temperature Tjv as shown by the deviation of 
the susceptiblity from linearity below the Weiss temper- 
ature. In the specific case of R=Y, oxygen stoichiomet- 
ric samples eventually order magnetically at low tem- 
perature (T7v=110K) due to structural distortion at T5 
> Tjv that lifts the three- fold degeneracy and Ocw—- 
508KQ. YBaCo407.o is thus a moderately GF magnet 



(f=6'(7n//Tjv 5). Interestingly, oxygen-rich samples re- 
tain three-fold symmetry and a disordered ground state 
down to the lowest temperature @. 

Although several neutron scattering studies have been 
recently conducted, the exact nature of the magnetic 
ground state in these systems remains unclear. Soda 
et al.^ have reported two magnetic transitions in single 
crystal of YBaCo407 both with short-range correlations, 
suggesting that the magnetic behavior of the triangular 
and Kagome layers are essentially decoupled. On the 
other hand, our powder diffraction experiments Q re- 
vealed that long-range magnetic order appears at a sin- 
gle transition (Tjv = 110K) for oxygen-stoichiometric sam- 
ple of YBaCo407. The magnetic structure derived from 
Rietveld refinement, shows that both magnetic Co sites 
bear a significant ordered moment, with a complex non- 
colinear magnetic configuration, contradicting the pic- 
ture of decoupled layers. 

Using single crystal neutron diffraction combined with 
Monte-Carlo (MC) simulations of the magnetic diffuse 
scattering collected at Qcw >T=130K> Tjv, we demon- 
strate that YBaCo407 represents a new GF system, char- 
acterized by the coexistence of long-range quasi-one di- 
mensional (ID) correlations and strong degeneracy in the 
two remaining directions. Our analysis reveals that the 
primitive magnetic clusters are corner-sharing triangular 
and bipyramidal units, and minimization of the intra- 
cluster exchange energy dictates configurations with zero 
total moment in each cluster. The form of the spin 
Hamiltonian uniquely promotes strong ID spin-spin cor- 
relations between apical Co(l) sites within chains run- 
ning along the c-axis. This contrasts with the fastly de- 
caying correlations along any direction in the a6-plane, 
characteristic of spin- liquid with 120° correlations. Our 
model provides an alternative explanation to that of 
Schweika et al. ^] who describe the magnetic scattering 
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from polycrystalline Yo.5Cao.5BaCo407 as arising from 
staggered chirality and for which the Co(l) site is non- 
magnetic. 

A single crystal of YBaCo407.o was grown by a float- 
ing zone technique in an optical image furnace. A densi- 
fied rod of the nominal composition was melted in a 20% 
02/Ar atmosphere, and the crystal grown at 1 mm/hr. 
Because of the high affinity for YBaCo407 to pick up oxy- 
gen on cooHng (ad, an after- heater was used to keep 
the growing crystal above 900°C. After growth was com- 
plete, the atmosphere was changed to 99.995 % Ar for 24 
hours and the crystal cooled down to room temperature 
over six hours. Oxygen content was verified by thermo- 
gravimetric analysis of a crushed sample of the crystal. 
Diffuse magnetic neutron scattering maps were recorded 
on the PRISMA time-of-fiight spectrometer at the ISIS 
facility, Rutherford Appleton Laboratory (UK). A ~ 200 
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FIG. 1: (Color online) Sketch of the magnetic topology in the 
YBaCo407 compound. Orange and blue spheres represent the 
inequivalent Co sites in, respectively, the Co(2) Kagome and 
Co(l) triangular layers (see text for details). The triangles in 
the Kagome planes are colored in blue for clarity (left panel) 
and bipyramids are colored in dark grey (right panel). Black 
and orange lines represent magnetic exchange interactions Ji 
and J 2 

mg piece of the single crystal of cylindrical shape (5 mm 
diameter, 2 mm height) was mounted on an aluminium 
pin with the [001] axis and [100] aligned vertically, in 
turn. For each axial orientation, scattering in the or- 
thogonal reciprocal-space planes were recorded by rotat- 
ing the crystal around the vertical axis. Data have been 
corrected for absorption and normalized to the incoherent 
scattering of a vanadium standard. The disordered mag- 
netic arrangements have been derived from direct Monte- 
Carlo (MC)_simulations [13] using a standard Metropolis 
algorithm [ll| using Heisenberg spins and isotropic ex- 
change interactions. Simulations have been performed 
using a supercell made of 12x12x12 conventional hexag- 
onal unit-cells and cyclic exchange boundary conditions, 
giving a total of 13824 spins. Magnetic neutron scat- 
tering intensities have been obtained by calculating the 



magnetic structure factor of the supercell [ig]- Powder- 
average was calculated by integrating the single-crystal 
scattering cross-sections in a hemisphere choosing dis- 
crete Q points that were homogeneously distributed. 

Neutron diffraction data collected at 330K are consis- 
tent with the trigonal space group P31c already proposed 
from measurements on powder samples Q. On cooling 
below Ts=300K, superlattice refiections corresponding 
to the orthorhombic Pbn2i symmetry 0, S], not shown, 
are observed for the three 120° domains. However for 
the sake of simplicity, we assume that the system can 
be modelled by a minimal set of exchange interactions 
on an ideal trigonal structure. There are four inequiva- 
lent Co-O-Co superexchange integrals to consider in this 
geometry, two in plane, between nearest neighbor Co(2) 
(non-capped triangle and capped triangle of the Kagome 
planes respectively), and two out of plane, between Co(2) 
and Co(l) (triangular layer above the Kagome layer or 
below) . A spin-dimer analysis using the Extended-Huckel 
Tight-Binding method ^] shows that the two in-plane 
and two out-of-plane interactions do not differ by more 
than 5%. We have therefore assumed only two inequiv- 
alent interactions Ji and J2, as shown in FiglTJ Fig[2] 
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FIG. 2: (Color online) Magnetic diffuse scattering for 
YBaCo407 in the hkO (left panel) and Okl (right panel) re- 
ciprocal planes. Each map shows the scattering cross sections 
obtained from Monte-Carlo simulations (left hand side) and 
the experimental magnetic scattering (right hand side) de- 
rived by subtracting the nuclear scattering collected at 330K 
from the total scattering at 130K. Dashed blue lines indicate 
the Brillouin zones (left) while solid lines (right) show the 
reciprocal lattice. 

shows the magnetic scattering in the hkO and Okl planes, 
obtained by subtracting the nuclear contribution (330K 
data) from the data collected at 130K, just above T^r. 
Although this procedure only partially removes the nu- 
clear scattering, as evidenced by the presence of residual 
sharp scattering at Bragg positions, it clearly highlights 
the magnetic signal. The latter is not resolution lim- 
ited and rapidly decays with the norm of the scattering 



vector (IQI), as expected from the niagnetic form fac- 
tor dependence. The in-plane magnetic scattering is rel- 
atively diffuse and concentrated along certain edges of 
the hexagonal Brillouin Zones (BZ). Remarkably, there 
is very little contribution in the first BZ and most of the 
intensity is found at the external perimeter of the second 
BZs. Reciprocal space cuts along directions perpendic- 
ular to the diffuse features leads to a full width at half 
maximum (FWHM) of O.lA'i, i.e. correlation lengths of 
about 30A. The nature of the diffuse scattering in the Okl 
plane is characterized by isolated diffuse spots at recipro- 
cal space positions 0,±i,2, the most intense features, and 
a spot a 0,|,0. The first set of peaks display anisotropic 
broadening with a FWHM along b* of 0.261-^ about 2 
times larger than the FWHM along c* of 0.12A~^, in- 
dicative of much longer correlations along the c-axis. 

MC simulations have been performed for a wide range 
of J2/J1 ratios The resulting phase diagram will 

be the object of a separate publication. It is found 
that the experimental magnetic scattering maps pre- 
sented here are best reproduced by a value of J2/J1 close 
to one. The corresponding calculated magnetic scatter- 
ing is plotted alongside the experimental data on Fig. 
[21 These results were obtained by simulated anneal- 
ing from a temperature T/J=100 down to a chosen fi- 
nal T/J. The maps in Fig[21 obtained by averaging the 
results of ten MC runs, are shown for T/J=4.0 which 
provides the best agreement with the experiment. This 
value is close to the calculated T/J=3.75 at T=130K if 
one take into account the mean-field AFM exchange in- 
teraction ((J) ~ 35K) extracted from the paramagnetic 
susceptibility The agreement with the data is ex- 
cellent, reproducing all main features, i.e. positions and 
relative intensities, in both scattering planes establish- 
ing that the essential physics is captured by our simple 
model. Crucial insights into the magnetic behavior of 
YBaCo407 are unveiled by parameters extracted from 
the zero-temperature ground state of the MC simula- 
tions. For example, analysis of the spin configurations, 
reveals that the magnetic lattice is made of triangular 
and bipyramid units, both with non-colinear arrange- 
ments and satisfying the a =X]i Si=0 rule. The latter 
point is illustrated in Fig[3] where the distribution of the 
total moment a taken over the entire lattice, shown for 
triangular and bipyramidal units, is centered around very 
small values. The small deviations from zero are related 
to inevitable residual fluctuations in the MC algorithm, 
frozen at low temperature. We also found that the mag- 
netic moments on the two apical positions of any bipyra- 
mid are systematically aligned parallel, a point developed 
in the next section. 

The CT=0 rule must be satisfied in any cluster of mutu- 
ally connected spins [l^ , such as for the triangular units 
imposing a angle of 120° between spins. However, spins 
in the trigonal bipyramid are not all mutually connected, 
since no nearest-neighbor exchange paths mediate inter- 
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FIG. 3: (Color online) Static spin-spin correlation functions 
for the Co(l) site along the z-axis (panel a) and in-plane 
([100], [010], or [l-10]-directions, panel b). Distributions of the 
magnetization within triangular units of the Kagome lattice 
(panel c) and within bipyramidal units (panel d). Sketches 
of the triangular and bipyramidal units together with typical 
spin configurations are depicted in the insets. 



actions between the apical sites labeled 4 and 5 in the 
inset of Fig (3] However, for a unique exchange integral 
J=Ji=J2, one can rewrite the Hamiltonian of the bipyra- 
mid up to a constant in a simple way: 



H = 



J 



^Si -2S4-Se 



1=1 



(1) 



This form also imposes (t=0 to minimize the energy. Ad- 
ditionally, Eq[l] indicates that only parallel alignment of 
S4 and S5 are possible ground states, in agreement with 
the configurations found in the MC simulations. This 
point is reinforced by examination of the static spin-spin 
correlation (SSC) function along the c-axis. The SSC for 
a given magnetic site Si in a direction R is defined as: 



(5(0)5(i?)) = ^ Si • Si+R 



(2) 



where the sum runs over all spins i in the lattice. In or- 
der to access spatial correlations at long distances in all 
directions, a MC simulation with a 24x24x12 supercell 
was performed. The SSC functions have been calculated 
for T/J=0.01 and are presented in Fig^ The SSC func- 
tions along the c-axis decay extremely slowly. This is 
illustrated for the apical Co(l) site, with 98% correla- 
tions at the tenth neighbor. This is consistent with Eq. 
1 since the bipyramids are corner-sharing along the c- 
axis and explains why the diffuse spots in the Okl map are 
sharper along c*. In contrast, correlations between Co(l) 
in-plane are lost much faster(25% at the ninth neigh- 
bor). They are also reminiscent of a 120°-configuration, 
characterized by the oscillating nature of the SSC with 
a periodicity R=3 and a value of -0.5 for R=l. Again, 
this behavior agrees with the topology of the lattice since 
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first neighbor bipyramids in the a6-plane are connected 
through a triangular unit that imposes such 120° config- 
uration. 

Further insight into the nature of the spin texture in 
YBaCo407 is provided by the correlations of both the 
vector (Kv) and scalar (Kg) chiralities for triangular 
units, as defined in [ij]- Indeed, in contrast to perfect 
Kagome lattices with staggered {q=y^3xy/3 ) or uniform 
(q=0) Kv, or to the non-zero Kg in the field- induced 
spin-canted structure seen in Kagome iron jarosite [l5| . 
our model doesn't exhibit any sign of such correlations. 

Finally, one can derive the powder-averaged magnetic 
cross section (Fig|4|) for several values of T/J includ- 
ing that which most closely matches our single crystal 
data. While the extracted intensity agrees with our neu- 



Co ions along the c-axis and short-range in the a6-plane, 
preserving 120° configuration, bear the signature of a new 
type of magnetic lattice. Despite the strong out-of-plane 
coupling between the Kagome layers in YBaCo407 acting 
as a field and constraining the spin configurations for half 
of the triangular plaquettes, the system retains a strong 
in-plane degeneracy. This unique feature should be rel- 
evant to explain physical properties for other analogs of 
the series and will no doubt spur experimental work in 
the samples where the three-fold symmetry remains at 
low temperatures. 

The authors would like to thank John Chalker for stim- 
ulating discussions. 
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FIG. 4: (Color online) Simulated magnetic scattering cross- 
section for a polycrystalline sample of YBaCo4 07 at vari- 
ous temperatures. The powder-averaged intensities have been 
calculated by numerical integration of the calculated single- 
crystal scattering cross sections (see text for details). 

tron powder data for YBaCo407 with the presence 
of a main peak at Q=1.35A~^, it also closely matches 
that observed in Yo.5Cao.5BaCo407 0], where the api- 
cal Co ions are reported to be non-magnetic. In that 
work, Schweika et al described the correlations in the 
latter system as arising from a q—y/3xy/3 structure, in 
which local zero-energy thermal weathervane defects, on 
the scale of one hexagon of the Kagome motif, results 
in an intense peak at 2qyy^=1.35A~^ but breaks longer 
range correlations, hence the weak signal they observed 
at q^y^—0.78A~^ . In our model, the peak at 0.78 
is naturally very weak near zero-temperature and is non 
existent at slightly higher ratios of T/ J due to the weak 
scattering within the first BZ. 

In summary, our magnetic neutron diffraction study in 
the spin-liquid regime of the frustrated YBaCo407 sys- 
tem, reveals that the magnetic lattice is best described 
as a 3D network of corner-sharing triangles and trigonal 
bipyramidal clusters, as illustrated on the right panel of 
FiglH Magnetic correlations, long-range between apical 
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